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ABSTRACT 

Context. The EROS-2 survey has produced a database of millions of time series from stars monitored for more than six years, allowing 
us to classify some of their sources into different variable star types. Among the so-called, long period variables (LPVs), which are 
known to follow particular sequences in the period-luminosity diagram, we find long secondary period variables whose variability 
origin remains a matter of debate. 

Aims. We analyze data for the 856 864 variable stars that have been discovered in the Large Magellanic Cloud (LMC) that are present 
in the EROS-2 database, to detect, classify, and characterize LPVs. 

Methods. Our method for identifying LPVs is based on the statistical Abbe test. It investigates the regularity of the light curve with 
respect to the survey duration in order to extract candidates with long-term variability. The period search is performed using Deeming, 
Lomb-Scargle, and generalized Lomb-Scargle methods, combined with a Fourier series fit. Color-magnitude, period-magnitude, and 
period-amplitude diagrams are used to characterize our candidates. 

Results. We present a catalog of 43 551 LPV candidates for the Large Magellanic Cloud. For each of them, we provide up to five 
periods, a mean magnitude in EROS-2, 2MASS, and Spitzer bands. Be - Re color. Re amplitude, and spectral type. We use infrared 
data to distinguish between RGB, O-rich, C-rich, and extreme AGB stars. Properties of our LPV candidates are investigated by 
ansilyzing period-luminosity and period-amplitude diagrams. 

Key words. Stars: AGB and post-AGB - late type- carbon - variables: general - Galaxies: Magellanic Clouds - Astronomical data 
bases: Catalogs 



1. Introduction 

Long period variables (LPVs) are red giant stars whose bright- 
ness varies on timescales from weeks to years. They are either 
mono- or multi -periodic, and include variable stars from the 
variability types formerly classified as Mira, semi-regular, and 
long secondary period Q^SP) variables. Among the distinctive 
characteristics of LPVs are the period-luminosity (PL) relations 
that they follow. A relation between the period and the absolute 
magnitude for Mira variables was first proposed by Gerasimovic 
(1928). Thanks to J, H, and K band near-infrared data, Glass & 
Evans (1981) were able to refine this relation that until then had 
displayed large scatter, as observations had been made using vi- 
sual magnitude. Wood & Sebo (1996) found a sequence in the 
(K, log P) plane that is parallel to the Mira PL relation at shorter 
periods and compatible with a first overtone pidsator population. 

The advent of large-scale surveys such as MACHO and 
OGLE has opened a new dimension in the study in general of 
variable stars, and in particular of LPVs. Initially devoted to 
the search of microlensing events, these surveys provided large 
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databases that made the statistical study of LPVs highly signif- 
icant. Wood et al. (1999) identified five sequences in the PL di- 
agram, using MACHO data (Alcock et al. 1992) of the Large 
Magellanic Cloud (LMC). They were labelled A, B, C, E, and 
D with increasing periods. The previously known PL relation 
of Mira stars falls on sequence C and corresponds to the radial 
fundamental mode of pulsation. Sequences A and B were iden- 
tified to be third, second, or first overtone pulsators and com- 
prise semi-regular variables. Ita et al. (2004), using OGLE-II 
(Udalski et al. 1997) and SIRIUS (Nagashima et al. 1999) sur- 
veys, showed that sequence B could be separated into two se- 
quences, called B and C. Using pulsation models from Wood 
& Sebo (1996) they concluded that sequence C is composed 
of "Mira-like" variables pulsating in their first overtone, com- 
plementing thereby the Mira fundamental mode pulsators of se- 
quence C. However, the theoretical periods for the second and 
third overtones do not fit their observed sequences B and A well. 
Therefore, models for these two sequences are still to be found. 

Theory predicts that fundamental radial mode pulsators have 
the longest periods, though sequences E and D appear to be at 
longer periods than sequence C. Since the study of Wood et al. 
(1999), sequence D is known to contain LSP variables. These 
stars were characterized as having two periods, with a period 
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ratio of the longer to the shorter period of 9 ±4, the longer period 
belonging to sequence D and the shorter, usually, to sequence B. 
Using OGLE-II and OGLE-III (Udalski et al. 2008) photometry, 
Soszynski et al. (2004b) suggested that sequence E may be the 
prolongation of sequence D at lower magnitude. As sequence E 
was found to be populated by ellipsoidal binaries, the shift of 
sequence E to periods shorter than those of sequence D would 
be due to the binary nature of the objects. The orbital periods 
of those stars should actually be twice as long as found, placing 
them on line with sequence D, which, as a consequence, would 
also contain binaries. However, Nicholls et al. (2010) excluded 
the hypothesis of binary systems for LSPs. Hence, the origin of 
variability for LSPs is still debated (Nicholls et al. 2009; Wood 
& Nicholls 2009; Nie et al. 2010). 

The analysis of the data from other large-scale surveys pro- 
vided further information about LPVs. While sequences A to D 
published by Wood et al. (1999) was assumed to consist of only 
asymptotic giant branch (AGB) stars. Kiss & Bedding (2003) 
demonstrated that the sequences extended below the tip of the 
red giant branch using OGLE-II and 2MASS (Skrutskie et al. 
1997) data, establishing the presence of red giant branch (RGB) 
pulsators as well. 

In this paper, we search for LPVs in the LMC and compile 
a catalog, from the database of variable stars provided by the 
EROS-2 survey. The EROS-2 survey was initiated in 1996, and 
during almost seven years monitored the Magellanic Clouds, the 
Galactic bulge, and four fields in the Galactic disk. Despite the 
rich database gathered during its seven years of operation called 
the EROS-2 catalog of variable stars that records, for example, 
856 864 vai-iables in the LMC fields (Mai-quette et al., 2011 in 
prep.), the survey remains to date under-exploited in the field of 
red-giant variable-star studies. 

The catalog published in this paper comprises 43 551 vari- 
able stars that we identify as LPV candidates in the EROS-2 
survey of the LMC. This survey is particularly well-suited as 
it covers a sky area of 88 deg^ on the LMC, compared to 13.5 
deg2 for MACHO and 40 deg^ for OGLE-III. 

We start Sect. 2 by recalling the main characteristics of the 
EROS-2 survey of the LMC and of its variable star content. 
The identification of LPV candidates of the LMC is then de- 
scribed in Sect. 3. The computation of pulsation periods and 
some characteristics of the selected LPV candidates are shown in 
Sect. 4. In Sect. 5, we discuss their infrared, period-luminosity, 
and period-amplitude properties. Conclusions are finally pre- 
sented in Sect. 6. 



2. The EROS-2 survey 

2.1. Presentation of the EROS-2 survey 

The EROS-2 survey operated from July 1996 to Mai'ch 2003 
with seven targets, namely the LMC, the Small Magellanic 
Cloud (SMC), the Galactic center, and fields around Mus, y 
Nor, y Set, and /3 Set in the disk of the Galaxy. Observations 
were made with the MARLY, aim Ritchey-Chretien telescope 
operating at La Silla Observatory in Chile. Two cameras, with 
a pixel size of 0.6" and a field of view of 0.7x1.4 deg^ on the 
sky, taken data simultaneously in the non-standard blue Be (420 
- 720 nm) and red Re (620 - 920 nm) bands. 

The observation properties of the EROS-2 LMC survey are 
summarized in Table 1 . Owing to some problems with the 'red' 
camera towards the end of the survey program, the time series 
in the R band are shorter than in the B band by about 300 days. 
Nevertheless, we use R band data for our analysis, since LPVs 
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Fig. 1. EROS-2 {Be — Re, Re) color-magnitude diagram of the ratio of 
the number of variable stars to the total number of stars. The color- 
magnitude diagram is divided into a matrix of 250x250 bins, each bin 
being 0.014x0.028 mag wide. The ratio is plotted according to a color- 
scale shown on the right of the figure. Only bins with at least ten stars 
and a ratio above 0.1 are plotted. Color version of the figures in this 
paper are available in the online version. 

are red and hence benefit from a higher photometric precision 
in the R band than in the B band. Data pre-reduction, dealing 
with corrections for a flat-field, bias, and dark cuiTent, was di- 
rectly made at the telescope. Data reduction and analysis were 
then made at the IN2P3 computing center in Lyon, France, with 
the dedicated programs package PEIDA (see Ansari 1996). Our 
work is based on the resulting light curves provided by the EROS 
team. 

A total of 87 million sources have been observed by EROS-2 
over its whole mission, 4.6 % of them being classified as variable 
(Marquette et al., 201 1 in prep.), which correspond to a database 
of time series for more than 4 million sources. One of the most 
populated target is the LMC, with a total of 28.8 million ob- 
served stars of which about 3 % were detected as variable (cf. 
Table 1). 

2.2. LMC variable stars in fhe EROS-2 database 

The ratio of the number of variable stars to the total number of 
stars across the color-magnitude diagram (CMD) is shown in 
Fig. 1. It displays the fraction of LMC stars if at least 10% of 
them are variable in each bin. Four structures are clearly visi- 
ble in this figure. The one around Be - Re - -0.25 mag and 
spreading from magnitudes Re = 18 mag down to Re = 14.7 
mag, represents the main sequence and would contain, among 
others, eclipsing binaries and variables of the Be type. The sec- 
ond one at Be - Re - 0.33 mag and Re - 18.8 mag is the RR 
Lyrae clump. Despite their faint luminosities, they are identifi- 
able thanks to their amplitude of variability of about magnitude. 
Another vertical structure is visible at Be - Re = 0.4 mag and 
from Re = 17 mag down to Re = 13.5 mag. This structure con- 
tains Cepheids, but might also be populated by some foreground 
stars from our Galaxy. The last structure is the AGB extending 
to the redder and brighter part relative to the tip of the RGB vis- 
ible at Be - Re - 11 mag and Re - 15 mag. Almost afl AGB 
stars are variable, and thus potential LPVs. The lack of points at 
the faintest magnitudes does not necessarily reflect an absence of 
variable stars, but rather an observational bias because the vari- 
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Fig. 2. Distribution of the time interval between two successive obser- 
vations of the star Im0010kl7507, in Re band, for which a total of 556 
good observations were recorded over 2061 days. 



ability of fainter stars is more difficult to detect at these magni- 
tudes owing to the larger photometric errors. 

2.3. EROS-2 time sampling of tine LMC observations 

The 88 fields constituting the EROS-2 LMC database were 
rather homogeneously observed during the visibility periods of 
the LMC throughout the survey. A typical distribution of the 
time interval between two successive observations is illustrated 
in Fig. 2 for the star Im0010kl7507, whose full light curve is 
shown in Fig. 4 (star nomenclature follows the rules defined by 
Derue et al. (2002)). The majority of observations are separated 
by two days and most of them by less than five days. The few 
long time intervals between 30 and 112 days represent gaps of 
from weeks to months in the data caused by the LMC not be- 
ing observable or without good enough air mass conditions from 
the EROS-2 site. This distribution of time intervals ensure that 
the EROS-2 survey is suitable for the detection of all kinds of 
known LPV stars, with typical periods ranging from 20 to 2 300 
days according to the GCVS (Samus et al. 2009). 

2.4. Data fiitering 

In principle, all data in the EROS-2 time series with a magnitude 
lower than 99.999 and an error value smaller than 9.999 should 
be assumed to be good. This typically results in about 389 good 
observations per star over 2 002 days in the Re band, and 477 
good observations per star over 2 324 days in the Be band. 

However, we apphed additional filters to the data in the fol- 
lowing way: 

- First, we kept only variables with Re <19 mag. These 

variables are expected to have reliable photometry, as the 
photometric precision reaches 10% at 19"" magnitude (cf 
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Fig. 3. NormaUzed distribution of 18609187 error values, on a log 
scale, from variables randomly taken, in Re band. 

Tisserand 2004). This allowed us to study variables as faint 
as the magnitude of the red giant clump. 

- Second, we removed all measurements performed during the 
first 50 days of the mission, since some of these data points 
have very small error values in the EROS-2 catalog, yet con- 
tain some obvious outUers. 

- Third, we eliminate in the remaining measurements the most 
extreme magnitudes around the mean magnitude by remov- 
ing one percent of the points from the highest and the lowest 
values. 

- Fourth, we eliminated all points with an error larger than one 
magnitude. We show however in the normalized distribution 
of more than 18 million error values from randomly taken 
variables (cf. Fig. 3), that only 0.64% of the points have an 
error larger than one magnitude. Removing points with large 
errors should not affect our analysis of LPVs, while it can 
remove outliers that may be present in the data. 

- Finally, we consider only sources with at least 40 points in 
their filtered Re hght curve. 



3. Identification of LPVs 

In general, it is found that LPVs have photometric variabiUty 
periods ranging from a few weeks to several years. Since the 
typical time interval between two successive EROS-2 observa- 
tions of the LMC is shorter than five days (see Sect. 2.3), the 
time series of LPVs are expected to vary smoothly with time 
in the EROS-2 database. We take advantage of this feature in 
the first step of our identification of LPV candidates by using 
a selection criterion based on the Abbe statistic (von Neumann 
1941) (Sect. 3.1). We then exclude in a second step non-LPV 
candidates from this list by using the color of the candidates 
(Sect. 3.1.3). 

3.1. The Abbe test 
3.1.1. Abbe statistic 

We use the Abbe statistic to identify smoothly varying time se- 
ries. We consider a time series (f;, x,), where x,- is the value at 
time ti of a measured quantity, with index ; running over the 
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Fig. 4. Time series of the LMC star Im0010kl7507 in the Rg magnitude. 
Times are in Heliocentric Julian day (HJD) - 2 450 000. 



number of measurements n. In our case, the measured quantity 
is the Re or Be magnitude. The Abbe statistic considers the dif- 
ferences (xi+i - X,) between two successive measurements, com- 
putes the sum of the square of those differences, and normalizes 
the result with the standard deviation of the time series and the 
number of observations. The Abbe statistic, which we hereafter 
simply call the Abbe value, is defined to be 



2(n-l) I.UiXi-x)^ 



(1) 



where x is the mean value of the measured quantity. The 
smoother the time series, the smaller the sum of the differences 
for a given standard deviation and number of observations. 

3.1 .2. Abbe statistic of irregularly sampled time series 

The value of Eq. (1) for a time series depends on the variability 
timescale of the light curve and the time sampling. In real data, 
the time sampling is not uniform. To assess the impact of the 
irregular time sampling on the Abbe statistic, simulations of time 
series of cosine functions = cos(2;r?,/P) were performed, 
with a time sampling typical of the EROS-2 observations of the 
LMC. We take for this purpose the time sampling of the star 
Im0010kl7507, whose light curve is shown in Fig. 4. 

Results are shown in Fig. 5. The simulations were performed 
with periods ranging from 0.1 to 4000 days. Each point in the 
figure represents the Abbe value for one periodic time series, 
where the period is plotted on the x-axis. From the top to bottom 
panel, an increasing level of noise is added to the time series. The 
error added to each point in the time series is computed from a 
gaussian distribution, of which the standard deviation o-e„ (noted 
as stdev(err)) relative to the peak-to-peak variability amplitude 
Avar of the error-free cosine values of the simulations (noted as 
ampl(rawTs)) is reported in each panel of the figure. We define 
the noise level as Ayar/o'err = ampl{rawT s) I stdev{err). When 
the noise level increases, the amplitude of the computed Abbe 
values becomes smaller. For the noisiest simulations, i.e. with 
'^var/o"err = 0.1, the time scries are so dominated by noise that no 
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Fig. 5. Values of the Abbe statistic for simulated cosine curves using 
the time sampling of the star Im0010kl7507 displayed in Fig. 4. The 
meaning of ampl(rawTs) and stdev(err) is given in Sect. 3.1.2 and that 
of stdev(ts) in Sect. 3.1.3. The x-abscissa is the period of the cosine 
in days, shown on a logarithmic scale. The vertical Une indicates the 
period at 20 days. 



correlation remains between the values of two successive points 
in the time series. In that case, the sum of the square differences 
(x,+i - x,)^ between two successive points equals, on average, 
twice the standard deviation, leading to a value of r ^ 1 .0. This 
conclusion is supported by the work of Strunov (2006), who 
shows that a normal distribution of uncorrelated values in a time 
series with more than 60 points leads to Abbe values that are 
normally distributed with a mean at 1.0 and with a variance of 



n-2 



' (K-l)(n+2)' 

The simulations with EROS-2 time sampling also demon- 
strate that the amplitudes of the oscillations in the Abbe values 
decrease towards shorter periods, which can be seen as a de- 
crease in dispersion around unity, and the Abbe values decrease 
towards zero for longer periods. This is important when selecting 
LPV candidates characterized by long periods. 



3.1 .3. Abbe statistic and long period variables 

From the conclusion of the section above, we can infer on a Umit 
to the Abbe value below which mainly variable stars with 
long variability timescales are selected. Fig. 5 shows that if 
is properly chosen, only a small fraction of the candidates with 
r < Him will have short periods, of either arotmd one day or, to 
a lesser extend, around 0.5 days (both due to the one day aUas- 
ing introduced in the time series obtained from a grotmd-based 
survey). 
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Fig. 6. Distribution of the ratio (r, J(rE,.r.Ri: on a log scale, for a sample 
20% of which are LMC EROS-2 variables that are randomly selected 
(172427 stars). 
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Fig. 7. Distribution of the Abbe statistic for a cleaned time series of all 
LMC variable stars (solid line). The dashed line shows the expected 
distribution of the Abbe values for simulated time series with normally 
distributed magnitudes. 



The Abbe value limit depends on the Avar/cerr level of the 
time series. Unfortunately, this quantity is unknown from obser- 
vations, since we do not know the variabiUty amplitude Avar of 
the noise-free light curve. Therefore, instead of using Avar/o'err, 
we use the ratio o-i^/o-f.„ where cti^ is the measurable standard de- 
viation of the observed time series (noted as stdev(ts) in Fig. 5). 
The mean value of this ratio over all simulated time series is 
reported in each panel of Fig. 5, and equals 7.1, 1.7, 1.2, and 
1.0 for Avar/ Cerr = 10, 2, 1, and 0.1, respectively. Fig. 6 shows 
the distribution of the ratio o-ts/o-f.„ for a sample of variables. 
The mean value of the distribution is equal to 5.89, hence our 
signal-to-noise ratios from the simulations are compatible with 
the values computed from the real data. The value of o-t^/osn for 
the star Im0010kl7507 is 9.92. If we consider the time sampling 
of this star to be representative of LMC LPVs in the EROS-2 
database, we expect all periodic stars with periods longer than 
20 days to have r < 0.4 (see top panel of Fig. 5). The population 
of stars satisfying this criterion would be only slightly polluted 
by periodic stars with periods of one day or its sub-multiples, as 
discussed previously. 

The distribution of Abbe values for all time series in the 
EROS-2 database of variable stars is shown as a soUd line in 
Fig. 7. Two populations are distinguishable in this distribution. 
The first one which peaks at r ^ 0.04, should correspond to the 



long timescale variable stars to which LPVs belong. The second 
population, around r = 1.0 corresponds to the time series that do 
not display a smooth variation in magnitude from one observa- 
tion to the next. The condition r < 0.4 suggested by the analysis 
of the Abbe statistic in Sect. 3.1.2 is seen in Fig. 7 to help us 
distinguish between the two populations to include as much of 
the first population as possible. 

The condition r < 0.4 on the Abbe value leads to the se- 
lection of 60723 candidates from the LMC variable stars, in 
which LPVs are to be found. We plot these candidates on the 
{Be - Re, Re) CMD of Fig. 8. We can safely assume that most of 
the red-giant variable stars matching our long-period variabiUty 
criteria were selected. On top of that, bluer variable types such as 
Be stars or long-period Cepheids for instance are also selected. 
We use the color of the red giant clump at Be -Re - 0.6 mag as 
the limit of the RGB, hence variable stars with a Be - Re > 0.6 
mag were chosen as potential LPV stars. The sample created 
consists of 43 583 EROS-2 sources. The criteria r > 0.4 and 
Be - Re > 0.6 mag selects 280 769 variable stars, indicating 
that variable stars redder than the red clump are not necessarily 
LPVs according to our criteria. A look at the hundred brightest 
light curves shows that they include: (1) binary stars, (2) stars 
with small amplitude and short variability, and (3) variables that 
have a long-term variability component associated with a sec- 
ondary variability component that is too short to result in a light 
curve smooth enough to give a small Abbe value. 



4. Period search 

4.1. Methods 

We search for periods characterizing the light curves of our 
LPV candidates using Deeming (Deeming 1975), Lomb-Scargle 
(Press et al. 2007), and non-weighted generalized Lomb-Scargle 
(Zechmeister & Kiirster 2009) methods. The Deeming method 
uses the Fourier transform applied to unevenly sampled light 
curves. The Lomb-Scargle method is derived from the least 
squares fitting of sine and cosine functions with a fixed mean 
value. Both methods are quite similar and should provide sim- 
ilar results, hence they are used as a double-check. The non- 
weighted generahzed Lomb-Scargle method adds a floating 
mean value computed from a sine function fitted to the time se- 
ries as an extra parameter to the classical Lomb-Scargle least 
squares. It prevents the selection of an incorrect period due to a 
poor phase coverage. For each of these period search methods, 
we use an iterative procedure to extract the significant periods 
that characterize a light curve. The standard deviation cr in the 
periodogram provides an estimate of the noise level. The signifi- 
cance of each peak in the periodogram is then evaluated relative 
to this standard deviation cr. A peak is considered significant if 
its value is at least 3 cr from the mean value of the periodogram. 
We follow a classical period-search procedure, i.e we first search 
for the most significant frequency in the periodogram of the orig- 
inal light curve (Fourier transform in the case of Deeming for ex- 
ample), fit the light curve with a sine function, and compute the 
residual from this fit. We then proceed similarly with the residual 
time series so obtained, and repeat the process until either no ad- 
ditional significant period is found or a maximum of twelve pe- 
riods are found. Figure 9 illustrates these steps with an example 
for the star Im0024n4304 in the R band. We show in the top pan- 
els periodograms resulting from the period search, the resulting 
model of the fit upon which we superimpose the raw light curves 
in the middle panels, and in the bottom panels the corresponding 
residual light curves after subtracting the fit from the raw light 
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Fig. 8. Number density in the EROS-2 LMC CMD for variable stars 
with an Abbe value <0.4. The vertical line at Be - Re = 0.6 mag shows 
the limit chosen for selection of LPV candidates. The CMD is divided 
into 250x357 bins of size 0.014x0.028 mag. The number of stars in 
each bin is plotted according to the color-scale shown on the right of 
the figure. 
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Fig. 10. Spectral window of the star lm0396114700. Top insert: zoom 
into the frequency interval used for period search. Bottom insert: 
Structure of the spectral window close to the one day alias. 



curve according to the period(s) found. We note that the periodic 
light curves that are not sinusoidal will have their harmonics ex- 
tracted during this process, so that the different periods may not 
be independent of each other. In our approach, those harmonics 
will be considered as separate periods. 

Numerically, the periodograms are computed for frequencies 
ranging from 1.818.10""^ to 0.1 day"', corresponding to periods 
from 10 to 5 500 days (2.5 times the mean EROS-2 duration) 
with frequency steps computed as 

Mt)-^, (2) 



^0.25 
p 0.2 

O 

10.15 

N 

§0.05 


2 4 6 8 10 

Number of periods 




where Ar is the time span of the time series and is an over- 
sampling factor that we set to 300. As most EROS-2 time series 
on the LMC have a duration around 2 000 days, we usually have 
adv(f)^ 2.10-''day-i. 

In Fig. 10, we show a typical spectral window for the LMC 
EROS-2 sources. Apart from the peak at the origin, the largest 
amplitude peaks are at frequencies of one day ' and higher. In 
the case of the period search for short-period variables (of pe- 
riod shorter than one day) these peaks will produce aliases in 
the power spectrum. In the bottom insert, we show the spectral 
window around one day '. We can see that the sampling allows 
us to clearly define the structure of the peaks. The two highest 
ones are at one solar day (1 day"'), and one sidereal day (1.0027 
day"'). As we applied a selection criterion to the Abbe values, 
only LPVs were retained. Spectral window for the frequency 
interval used for the period search for our LPV candidates is 
shown in the top insert of Fig. 10. The peak clearly visible at a 
frequency of 0.0027 day"' is produced by the one-year artifact. 
Its normalized amplitude is smaller than 0.25. The aliases cre- 
ated in the periodogram then, are not expected to be significant. 
According to the values of the peaks found at frequencies lower 
than 0.1 day"', we can assume that the period(s) found for our 
LPV candidates are not produced by aliases. 



Fig. 11. Normalized histogram of the number of periods found for the 
Deeming, Lomb-Scargle and generalized Lomb-Scargle methods. 



4.2. Results 

Figure 11 shows a normalized histogram of the number of 
periods found for each of the three methods, according to the 
criteria discussed in the previous section. Each star appears 
only once, which means for example that a star for which three 
periods are found is not counted in the two-period bin. The 
proportions of variable stars with the same number of periods 
appear to be comparable independently of the method used. 
However, for roughly six percent of stars no period was found 
with the Lomb-Scargle and the generalized Lomb-Scargle 
methods, whereas Deeming found at least one. We select these 
1 244 variables stars and show the histogram of their first period 
found with the Deeming method in Fig. 12. We can see that 
for the majority of these stars the period exceeds 2000 days, 
which is the mean EROS-2 time length. Thus, a complete phase 
coverage is impossible in these cases. Nevertheless, visual 
inspection of these stars shows that most of them exhibit clear 
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Fig. 9. Period search example for star Im0024n4304. The top panels show the successive power spectra of the raw light curve (left panel) and of 
the first residual (right panel). The middle panels display the light curve in R band (black points with error bars), with the fitted Fourier series 
superimposed (red curve) for the first (left panel) and first and second (right panel) periods found. The bottom panels show the residuals after 
subtraction of the fits from the light curve. 
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Fig. 12. Histogram of the first period found with the Deeming method, 
for variables for which a period was found with neither the Lomb- 
Scargle nor generalized Lomb-Scargle methods. 



long-term variability features. 

For each period-search method, we consider up to five pe- 
riods. If three period values, each found by a different method, 
are within a range of +10% we define a new period value that is 



equal to the mean of the three corresponding periods determined 
by the different methods. 

For almost 14% of our LPV candidates, no common period 
was found by all three methods. Visual inspection of these time 
series shows that for about 10% of them we find a typical LPV 
light curve. A cause of their exclusion from the sample is that the 
highest peak in their periodogram is at an harmonic period for 
one of the methods, while in other cases it is at the fundamental 
period (3% of the cases). We reclassify these stars by choosing 
the period of the fundamental mode. Figure 13 shows the peak- 
to-peak amplitude in the band versus the standard deviation 
in the error values for the remaining 90% of the stars where no 
clear LPV feature was observed. Most of them appear to have a 
small standard deviation in the error values (o-grr < 0.05), even 
if a small group of variables is visible at large standard deviation 
with a larger mean amplitude (crg,.^ > 0.1 and Re Amplitude> 
1.5 mag). A visual inspection of the 72 variables with o-g,.,. > 
0. 1 reveals that 47 of these stars are affected by data reduction 
problems. As a consequence, many outliers are visible in these 
light curves, which explains the higher amplitude values. These 
outliers have errors smaller than expected for their magnitude, 
which cause the larger standard deviation value in the errors. 
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Fig. 13. Peak-to-peak amplitude versus the standard deviation in the 
error values. Stars without any period simultaneously found with the 
three-period search methods, are shown as black points. Red crosses are 
stars showed in Fig. 12. Blue triangles are RGB and DY PER identified 
by Tisserand et al. (2009). 



Table 2. Results of cross-identifications between EROS-2, OGLE III, 
and MAGHO LPV catalogs. 



Catalog 


Number of 


Common LPVs 


LPVs found 


LPVs found 




LPVs in cat. 


(EROS-2, cat.) 


in EROS-2 only 


in cat. only 


OGLE III 


91995 


27 940 


15 643 


64054 


MACHO 


56453 


25 812 


17771 


30641 



The resulting variability is then very unreliable. We therefore 
exclude them from further study. 

Visual inspection reveals that some non-periodic variables 
are present in our sample. We identified 32 variables of R 
Coronae Borealis (RCB) and DY Per types already known in 
EROS-2 and listed in Tisserand et al. (2009). We show them in 
Fig. 13 as triangles. These stars are flagged in the online cat- 
alog containing 43 583 EROS-2 sources (see Sect. 5.2) but are 
not taken into account for the following study. In Fig. 13, we 
also show, as red crosses, the 1 244 previously identified vari- 
ables for which only the Deeming method find period(s). These 
stars lie in the group with the smaller error dispersion, which 
lends support to the assumption that variables with small (Terr 
are probably genuine variables with a very long period. 
We then obtained 6 008 variables without finding any period 
and 37543 variables (86% of the LPV candidates) for which at 
least one period value is retrieved, distiibuted as follows: 9666 
variables with a single period retrieved, 9 904 with two periods, 
8 595 with three periods, 5 845 with four peiiods, and 3 533 with 
five periods. Finding LPVs with such a large number of periods 
is unexpected, but is probably caused by our keeping harmonics 
as independent periods. We created a catalog with 43 551 LPV 
candidates (see Sect. 5.2). The sample of LPVs studied in the 
following sections is composed of the 37 543 variable stars for 
which at least one period was found. 

We performed cross-identifications of all our LPV candi- 
dates with two LPV catalogs available from both the OGLE III 
(Soszynski et al. 2009) and MACHO (Eraser et al. 2008) sur- 
veys. As some of the EROS-2 RCB and DY Per stars are consid- 
ered as LPVs in both the OGLE-III and MACHO catalogs, we 
kept them for these cross-identifications. A two arcsec search 
radius was used to identify counterparts to our LPV candidates 
in these two catalogs. For cases where one LPV candidate gives 
two matches, we keep the one with the closest coordinates. The 
results of this cross-identification are summarized in Table 2. 



We derived 64% of our LPVs from the OGLE III catalog. 
Most of the remaining 36% of LPVs that were found only in 
EROS-2 are in fields of the LMC not covered by OGLE III. 
The large number of LPVs found only in OGLE III correspond 
mainly to OGLE small-amplitude red-giant variables (OS ARGs) 
(Wray et al. 2004), which constitute 86% of the OGLE III LPV 
catalog. As shown in Soszynski et al. (2004a) (see their Fig. 2), 
OSARGs can have amplitudes as small as few millimagnitudes 
in I band. Since Fig. 3 shows that most of EROS-2 errors have 
values between 0.01 and 0.025 mag in Re band, an identification 
of all OSARGs as vaiiable stars in EROS-2 is not very likely. 

We found that 59% of the LPVs are in common be- 
tween EROS-2 and MACHO, this smaller percentage com- 
pared to OGLE III being caused by a smaller sky coverage 
for the MACHO survey. Some of the 30641 LPVs found 
only in MACHO may come from small amplitude LPVs, since 
6 031 stars have an amplitude in R band smaller than 0.025 
mag. Another part would consists of the 11215 LPVs from 
the MACHO catalog affected by the one-year artifact (cf. 
Eraser et al. (2008)) that are not selected with EROS-2 (See 
Sect. 4.3.1.). The sky coverage is an important factor as well. 
Most of the LPVs found only in EROS-2 are in fields of the 
LMC not covered by MACHO, mainly external fields but also 
those in the LMC bar. 

4.3. Characterization of LPVs 

4.3.1. Mono-periodic LPVs 

More than 22% of our LPV candidates are found to have only 
one significant period according to the criteria mentioned in the 
two previous sections (i.e a period with a significant peak value 
in the periodogram and found to be in agreement with the three 
peiiod search methods). The Abbe values are plotted against the 
peiiod with color-coded amplitude in Fig. 14. We first note that 
no variables have a small Abbe value (less than 0.1) and short 
period (shorter than 50 days), as expected from Sect. 3.1.2. 

The distribution of periods in the top panel of Fig. 15 shows 
an increase in the number of mono-periodic LPVs up to log 
P = 2.3 (P=199.5 days). Then, instead of a monotonic de- 
crease of the distribution towards higher periods, two peaks are 
visible at log P - 2.51 and log P = 2.6 (P=323.6 days and 
P=398.1 days). Raw EROS-2 time series contain a one year ar- 
tifact caused by the seasonal effect, which is visible as linear 
trends in magnitude. In the light curves provided, this artifact 
was corrected to some extent (cf Tisserand 2004) and explains 
the smaller number of stars with a period of around one year, as 
indicated by the vertical dashed line in Fig. 15. Variables with 
peiiods a few weeks shorter or longer than this one year artifact 
could still contain these trends in their light curves, hence pro- 
duce an artificially larger number of LPV candidates with peri- 
ods around one year. 

4.3.2. Double-periodic LPVs 

For double-periodic variables, the middle and bottom panels of 
Fig. 14 display the same general characteristic as mono-periodic 
LPVs. The distribution of periods in the bottom panel of Fig. 15 
also shows a larger number of stars close to the one-year artifact, 
especially for the first peiiods slightly above one year visible as 
a strong peak at log P - 2.62. The distribution of secondary 
periods is more uniform, mainly for periods between 40 days 
(log P - 1.6) and one year. 
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Fig. 14. Top panel: Abbe values versus period on a log scale, for mono- 
periodic LPV, with color-coded amplitude. Middle panel: Abbe values 
versus first period on a log scale, for double-periodic LPV, with color- 
coded amplitude. Bottom panel: Abbe values versus second period on a 
log scale, for double-periodic LPV, with color-coded amplitude. 

5. Main properties of the sample 

As LPVs are known to follow distinctive sequences in the 
period-luminosity diagram, we identify and discuss here the 
properties of these sequences for our LPV sample. The surface 
abundances of red giant stars can be dominated by either oxy- 



300 




log(pcriod) (days) 



Fig. 15. Top panel: Histogram of the periods for mono-periodic LPVs, 
on a log scale. Bottom panel: Histogram of the periods for double- 
periodic LPVs, on a log scale. The black line is for the first periods 
found, the red one is for the second periods found. Vertical dashed lines 
at log P = 2.56 indicate the period at one year. 



gen (O-rich) or carbon (C-rich), depending on their evolution- 
ary phase. Some of them are obscured by circumstellar dust. 
Visual photometry, then, does not allow a direct comparison of 
these stars with the lightly obscured ones. Hence, near- and mid- 
infrared observations are needed to investigate the properties of 
RGB and AGB stars. Thus, in addition to EROS-2 data, we use 
2MASS and Spitzer infrared data for the analysis of our candi- 
dates. 

5.1. Infrared properties 

We cross-matched the 43 551 EROS-2 LPV candidates with the 
2MASS All-Sky Catalog of Point Sources (Cutri et al. 2003) us- 
ing the VizieR Services at "Centre de Donnees astronomiques 
de Strasbourg". Using a search radius of 1 arcsec, we matched 
41 502 stars (95.3% of the LPV candidates) in J and K, bands, 
with a mean difference between EROS-2 and 2MASS coordi- 
nates of 0.31 arcsec, and a cr,/=0. 16 arcsec. For three of these 
LPV candidates, the cross-matching ends with two possible 
2MASS sources. In these cases, we kept only the source with 
the coordinates closest to those of our candidate and with the 
best 2MASS quality flag. 

We also cross-matched our candidates with the Spitzer IRAC 
data at 3.6 microns, hereafter [3.6], from the SAGE database 
(Meixner et al. 2006). Within a search radius of 1 arcsec, 36 847 
stars (84.6% of the LPV candidates) were cross-matched, with a 
mean distance between EROS-2 and IRAC coordinates of 0.33 
arcsec, and a crj=0.18 arcsec. For three of these LPV candidates, 
we obtained two possible IRAC sources. In these cases, we kept 
only the source with the coordinates closest to those of our can- 
didate. 

Several methods exist to help us distinguish between O-rich 
and C-rich stars. Among these are the use of the reddening-free 
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Fig. 16. Number density of the EROS-2 LPV candidates cross-matched 
with the 2MASS catalog, in the J-A:.,, K, CMD. The CMD is divided 
into 200x200 bins with a bin size of 0.061x0.043 mag. 
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Fig. 17. Histogram of J-Ks for variables with 10 < Ks (mag)< 1 1 from 
Fig. 16. 



Wesenheit index W/ combined with periods (Soszynski et al. 
2005), the use of either {J - K,V - I) color-color diagram, the 
(WjK, Wj) diagram (Soszynski et al. 2009), or the variability 
slope parameter a,, versus color amplitude ampvi (Wisniewski 
et al. 201 1). Another way to distinguish between different spec- 
tral types in red giants, is to look at the J-K color index, O- 
rich stars usually having a lower (J - K) value than C-rich (cf. 
Nikolaev & Weinberg 2000). Fig. 16 shows the {J-K,, K,) Hess 
diagram for the LPV candidates. The two highest density zones 
are clearly separate at Ks=l2 mag, which is the magnitude of 
the tip of the red giant branch (TRGB). We defined four zones in 
that diagram as follows: 

- Zone 1 is mainly composed of RGB stars, although a con- 
tamination by faint AGE stars close to the TRGB is possible, 
as shown by Cioni et al. (2006). 

- For Zone 2, we use Cioni et al. (2006) equations [6] and [7] 
to separate O-rich from C-rich AGE stars, respectively zones 
2 and 3 in Fig. 16, assuming a metallicity of Z=0.008 for the 
LMC. 
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Fig. 18. Number density of the EROS-2 LPV candidates cross-matched 
with the IRAC catalog, in the J- [3.6], [3.6] CMD. The CMD is divided 
in 250x250 bins with a bin size of 0.04x0.037 mag. 
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Fig. 19. Histogram of 7 - [3.6] for variables with 9 < [3.6] (mag)< 10 
from Fig. 18. 
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Fig. 20. Histogram of [3.6] for variables with 1 < J - [3.6] (mag)< 1.5 
from Fig. 18. Dashed line at [3.6] = 11.85 indicates the magnitude of 
the TRGB. 



Zone 2, according to their models, contains O-rich thermally 
pulsating AGB (TP -AGB) stars, with presumably the bright- 
est part (K, < 10.6 mag) populated by young M supergiants 
and luminous O-rich M stars as discussed in Blum et al. 
(2006). 
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Fig. 21. (Be - Re, Re) CMD of the EROS-2 LPV candidates, according 
to the classification scheme into RGB, 0-rich, C-rich, and extreme AGB 
stars defined in Sect. 5. 



- Zone 3 contains C-rich TP-AGB stars, which evolved from 
the 0-rich, though a mixture of these two populations is pos- 
sible either by C-rich stars in zone 2 or by O-rich stars in 
zone 3. Fig. 17 shows the histogram of J - for variables 
from Fig. 16 with 10 < (mag)< 11. According to our def- 
initions of the zones, the highest peak is due to O-rich AGB 
stars from zone 2. For J - > 1.4 mag, C-rich AGB stars 
form the predominant population up to J - = 2.1 mag. 

- Zone 4 is composed of stars with J - > 2.1 mag, which 
tend to become fainter as they become redder, in constrast 
to O-rich and C-rich stars from zones 2 and 3. This popula- 
tion of stars called extreme AGB stars are known to be dusty 
evolved AGB stars (see Vijh et al. 2009). 

Owing to the presence of circumsteUar dust from extreme 
AGB stars, we require mid-infrared data to allow the compar- 
ison of this population with less obscured O-rich and C-rich 
AGB stars. This is illustrated in Blum et al. (2006) and Ita & 
Matsunaga (2011), where the use of Spitzer SAGE data places 
extreme AGB stars on the brighter end of the C-rich branch, 
brighter stars being the redder ones. Combining 2MASS and 
SAGE data, Blum et al. (2006) showed that the (/ - |3.6|, |3.6|) 
CMD is the most convenient to differentiate between O-rich, C- 
rich, and extreme AGB stars. In Fig. 18, we plot the 36 220 LPVs 
(83.1% of the LPV candidates) for which we have both J and 
[3.6] data, as a (/ - [3.6], [3.6]) Hess diagram. We find again 
the same features as in Fig. 16. As expected, extreme AGB stars 
are generally brighter than C-rich AGB stars. Furthermore, C- 
rich stars are brighter than O-rich ones. This is clearly visible in 
the histogram in Fig. 19, showing the distribution of 7 - [3.6] 
for variables from Fig. 18 with 9 < [3.6] (mag)< 10, where 
C-rich stars are the main population, and the remaining O-rich 
stars are M supergiants. The signature of the TRGB is visi- 
ble on the distribution of the [3.6] magnitude for LPVs with 
1 < / - [3.6] (mag)< 1.5 in Fig. 20 as a minimum of the dis- 
tribution at [3.6] - 11.85 mag. As no models are available for 



(7 - [3.6], [3.6]) diagram, we keep Cioni's models to select O- 
rich and C-rich TP-AGB stars from Fig. 16, while extreme AGB 
stars are selected according to their /- [3.6] color and RGB stars 
with respect to their [3.6] magnitude. For the following study, 
LPV candidates are thus classified as follows: 



- RGB stars: LPV stars that are fainter than the TRGB, i.e 
[3.6] > 11.85 mag. There are 17427 stars of this type in 
our sample. 

- Extreme AGB stars: those with /-[3.6] > 3.1 mag, cf Fig. 19 
and brighter than the TRGB in [3.6] band. These are the less 
numerous stars in our sample with 859 stars found. 

- O-rich stars belong to zone 2 from Fig. 16, and not already 
identified as RGB or extreme AGB stars. These are the most 
numerous AGB stars in our sample with 10 807 stars found. 

- C-rich stars belong to zone 3 from Fig. 16, and not already 
identified as RGB or extreme AGB star. There are 6612 
LPVs of this type in our sample. 



Figure 21 shows the distributions of these four groups of stars in 
the (Be -Re, Re) CMDs. We can see that it would be much more 
difficult to separate the four groups with EROS-2 data alone. 
The TRGB for example is highly contaminated by O-rich and 
C-rich AGB stars in the (Be - Re, Re) CMD. It is also impossi- 
ble to disentangle O-rich from C-rich stars as they have the same 
Be - Re color and Re magnitude range. Extreme AGB stars ap- 
pear fainter than other stars, as expected for stars with thick cir- 
cumsteUar dust shells, but not redder. This diagram confirms the 
interest of using infrared data to distinguish LPVs with different 
abundances. 



5.2. Catalog 

In our catalog, we present the 43 583 EROS-2 sources extracted 
as LPV candidates based on the Abbe value criterion and red 
clump color (cf. Sect. 3.1.3). Table 3 shows a sample of the cat- 
alog and contains the following information: 
Column 1: EROS-2 ID of the star according to the definitions by 
Derue et al. (2002). 

Columns 2-3: equatorial coordinates, RA and DEC, for the 
epoch J2000. 

Column 4: mean Re magnitude. 

Column 5: mean Be - Re color. 

Column 6: Re band amplitude as a peak-to-peak value. 

Columns 7-11: up to five periods, in days, for the Re band, if 

found with all the three period search methods. In the following 

analysis we will refer to these periods. 

Column 12: the number of periods found and given in columns 
7-11. The value is equal to in case no common period was 
found by all the three period search methods. 
Column 13: 2MASS J magnitude. 
Column 14: 2MASS magnitude. 
Column 15: Spitzer [3.6] magnitude. 
Column 16: the classification as defined in Sect 5.1. 
Column 17 :& flag value. Non-periodic RCB or DY Per type stars 
from Tisserand et al. (2009) are flagged with value 91 or 92, 
respectively. Stars found to have harmonic periods among se- 
quences D and E are flagged with the value 1. Flag is set to in 
all the other cases. 

The full catalog is available on-line at the CDS. 
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Table 3. Data sample of the LPV in the LMC from EROS-2. See Sect. 5.2 for column description. 



ID 

EROS-2 


RA 
(J2000) 


DEC 
(J2000) 


(mag) 


(Bi -«£),„„„ 
(mag) 


Amp.RE 
(mag) 


(days) 


P2r, 
(days) 


P^R, 

(days) 


(days) 


P5r, 
(days) 


NPer 


J 

(mag) 


(mag) 


[3.61 
(mag) 


Spectral 
type 


Flag 


Im0020n24738 


82.63144 


-69.35374 


13.923 


1.232 


0.333 


586.4 














1 


12.096 


10.7 


10.24 


C-rich 





Im0020n25074 


82.53326 


-69.35603 


16.832 


0.647 


0.258 




















999.9 


999.9 


15.213 


N/A 





Im0020n25567 


82.59598 


-69.35854 


14.826 


1.082 


0.124 


30.2 














1 


13.33 


12.241 


12.14 


RGB 





Im0020n25716 


82.26153 


-69.36008 


14.762 


1.303 


0.218 


46.8 


404.6 


2782.0 


35.5 


50.8 


5 


13.255 


12.116 


11.893 


RGB 





Im0020n26052 


82.27126 


-69.36201 


14.032 


1.550 


0.514 


625.0 


326.0 


88.3 


83.7 





4 
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Fig. 22. (log P, ATj) diagram for all the periods of LPV candidates. The Fig. 23. Analogous to Fig. 22 for the [3.6] magnitude versus log of the 
black points are RGB stars, blue points are 0-rich AGB stars, green periods, 
points are C-rich AGB stars, and red points, extreme AGB stars. 



5.3. Period-magnitude relations 
5.3.1. Identification of sequences 

Figure 22 shows the (logP, Ks) diagram of our LPVs for all the 
periods, which means for example, that a variable with three pe- 
riods will be represented by three points in the diagram. We color 
code the diagram according to the RGB, O-rich, C-rich, and ex- 
treme AGB star classification (see Sect. 5.1). Six sequences are 
visible. Five of them, namely A, B, C, D, and E, have been stud- 
ied in the seminal paper of Wood et al. (1999). Ita et al. (2004) 
subdivide the sequence B into the B and C sequences, using 
OGLE-II and SIRIUS data. We follow the same notation for our 
sequences as in these two works. 

Extreme AGB stars clearly produce a large scatter and a mix 
of population in the sequences. As shown in Riebel et al. (2010) 
and Ita & Matsunaga (201 1), the scatter in the period-magnitude 
relations is smaller when using mid-infrared data. Comparing 
2MASS and SAGE data, Riebel et al. (2010) showed in (log P, 
magnitude) diagrams that the sequences are more clearly defined 
when using the Spitzer [3.6] band. The (log P, [3.6]) diagram in 
Fig. 23 indeed, displays extreme AGB stars at the bright end, 
aligned with C-rich star sequences. Globally, the sequences are 
more easily differentiated using the [3.6] band than the band. 
Hence, we use the (log P, [3.6]) diagram for further analysis. 



In the (log P, [3.6]) diagram in Fig. 24, we show all the pe- 
riods found for our LPV candidates (including stars without de- 
fined spectral type). The signature of the TRGB is easily visible 
as a drop in star density at [3.6] = 11.85 mag for the entire pe- 
riod range, as seen in other works (cf. Kiss & Bedding 2003). 

We defined sequence boundaries by firstly connecting points 
in the valley of lowest density between the sequences. These 
points were identified as local minima of the distributions of 
log P for three magnitude intervals: 12 < [3.6](mag)< 12.5, 
11 < [3.6](mag)< 11.5, and 10 < [3.6](mag)< 10.5. Secondly, 
as the density of points is too low at the bright and faint mag- 
nitudes we visually add points at [3.6] = 8 mag and [3.6] = 14 
mag to adjust the shape of the boundaries to the sequences. We 
performed two least squares fits for each of these six sequences, 
one for the region of the diagram brighter and one for the re- 
gion fainter than the TRGB. The parameters of the equations we 
obtained are summarized in Table 4. Except for sequence A, all 
the sequences show a steeper slope for AGB stars than for RGB 
ones. 

5.3.2. Cleaning of harmonics 

In the following sections, we analyze mono and double periodic 
LPVs. Since light curves of ellipsoidal variables have minima of 
alternating depth, our period search could lead to the selection 
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Fig. 24. (log P, [3.6]) diagram for all the periods selected for the LPV 
candidates. The dashed line indicates the [3.6] magnitude of the TRGB. 
The red lines delimit the sequences as described in Sect 5.3.1. 
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Fig. 25. (log P, [3.6]) diagram for LPVs with only one period identified. 
The black points are RGB stars, blue points are 0-rich AGB stars, green 
points are C-rich AGB stars, and red points, extreme AGB stars. 




Table 4. Linear fit parameters for the sequences from Fig. 24 with re- 
spect to the TRGB. 
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RGB 


AGB 




Slope 


Intercept 


Slope 


Intercept 
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-5.64 


20.31 


-4.00 


17.77 
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-3.52 


17.93 


-4.70 


19.76 
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18.95 
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20.64 
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-2.41 
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-4.33 


20.72 


E 


-2.77 


18.72 


-3.61 


20.31 


D 


-3.05 


20.36 


-3.35 


20.71 



of harmonic periods. We handle these harmonics for double- 
periodic LPVs that are potential binaries, as follows: if a LPV 
has one of its periods in sequence D and the other in sequence 
E, with Pd - 2Pe+IQ% we remove the sequence D period and 
attribute only the sequence E period to the LPV. This means that 
we have binaries with minima of alternating depth and binaries 
with minima of similar depth in sequence E. A second step was 
to study all the LPVs that remain in the sequence D. We looked 
by eye at the raw and phase-folded light curves. If it appears to be 
typical of an ellipsoidal variable, we attribute a period P - Pq/I 
as their only period. For LPVs with more than two periods found, 
we attribute a flag in the catalog if two of these periods are mul- 
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Fig. 26. Cumulative distribution for LPVs with only one period identi- 
fied for the sequences in Fig. 25. Each line represents a sequence, which 
are, from the most populated (top line at [3.6] = 16.5 mag) to the least 
populated (bottom line at [3.6] = 16.5 mag): Seq. E, Seq. C, Seq. C, 
Seq. D, Seq. B, Seq. A. 



tiple by a factor of two, with one period in sequence D and any 
other one in sequence E. 
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5.3.3. Mono-periodic LPVs 

Data in the [3.6] band were available for 7 532 out of 9 666 LPVs 
classified as mono-periodic. Fig. 25 displays the (log P, [3.6]) 
diagram for these stars. The most populated sequence for the 
mono-periodic LPVs appears to be sequence E, as shown in the 
cumulative distribution of Fig. 26. Sequence E is known to con- 
tain ellipsoidal binaries, effectively found with only one period, 
but with half the value of the orbital period because of the meth- 
ods used for the period search, as shown in studies of binary sys- 
tems (see Soszynski et al. 2004b; Derekas et al. 2006). As visible 
in Fig. 25, variables from this sequence are present at all mag- 
nitudes. Nevertheless, Fig. 26 shows that most of the sequence 
E variables are fainter than the TRGB ([3.6] = 11.85 mag) in 
contrast to all the other sequences. The second most populated 
sequence, with more than 25% of the mono-periodic LPV sam- 
ple, is the sequence C. The stars on this sequence are radially 
pulsating in the fundamental mode, the so-called Miras. Except 
for a very few cases in sequences C and D, sequence C is the 
only one containing extreme AGB stars. 

Some variables with only one period are distributed over 
sequences C, B, and A, populated by overtone pulsators. 
Sequence C is populated by first overtone pulsators (cf. Wood 
& Sebo (1996), Ita et al. (2004)). Sequences A and B are ex- 
pected to be populated with semi-regular small amplitude red- 
giant variables pulsating in their second or third overtones. The 
nature of these variables is not easy to determine, because avail- 
able models of low-mass (in the range 0.8 to 1.5 Mq) pulsating 
stars in overtones modes are unable to reproduce the sequences 
observed (cf. Wood et al. 1999). These sequences are proba- 
bly composed of a mix of classical semi-regular variables and 
OSARGs that may require different physical processes to ex- 
plain their variability (Soszynski et al. 2007). 

Only a few stars are visible in sequence D, which usually 
contains LSP variables, that have more than one period by defini- 
tion. Visual inspection of our LPV candidate light curves in this 
sequence reveals that it contains mostly stars that exhibit a LSP 
feature with a second period that is too short (i.e. shorter than 
10 days) to be detected by our period search. Some other types 
such as in-egular variables with a long-term trend, and peculiar 
stars exhibiting long period variability with an amplitude varia- 
tion, probably due to spotted stars (cf. Marquette et al. (2008)), 
are also found in this sequence. 



5.3.4. Double-periodic LPVs 

The most numerous subtype in our sample are double-periodic 
LPVs, making up one quarter of LPV candidates (cf. Table 5). 
For each of these stars, we assign two sequences with respect to 
the position in the (log P, [3.6]) diagram of the periods found. 
Figure 27 shows the (log P, [3.6]) diagram of the first period 
found, for the double-periodic LPVs. One can see the few ex- 
treme AGB stars in sequence C relative to the mono-periodic 
case (cf. Fig. 25). Figure 28 (top panel) shows the cumulative 
distribution for each sequence in the [3.6] band. It shows that 
the majority of the population consists of LPVs in sequence D, 
supposedly LSPs, for which almost half of these stars are fainter 
than the TRGB. Double-periodic stars with their first period 
in fundamental (sequence C) or low-order overtone mode (se- 
quence C) are the next most numerous variables. There is still a 
noticeable population of stars in sequence E. Semi-regular vari- 
ables from sequences A and B are more numerous than for the 
mono-periodic case (+44% of LPVs for sequence A and -1-69% 
for sequence B), but remain the least populated sequences. 
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Fig. 27. Log of the first period versus [3.6] diagram for LPVs with two 
identified periods. Colors are as in Fig. 22. 
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Fig. 28. Cumulative distribution for double-periodic LPVs for the se- 
quences in Fig. 25. Each line represents a sequence, which are, from 
the most populated (top line at [3.6] = 16.5mag) to the less populated 
(bottom line at [3.6]=16.5mag): For the first period (top panel): Seq. D, 
Seq. C, Seq. C, Seq. E, Seq. B, Seq. A. For the second period (bottom 
panel): Seq. C, Seq. D, Seq. B, Seq. C, Seq. E, Seq. A. 



Figure 29 shows the (log P, [3.6]) diagram of all the second 
periods found. In Fig. 30, we show six (log P, [3.6]) panels for 
the second period according to where the first period was found. 
This allows us to see the origins of the sequences of the (log 
P, [3.6]) diagram in Fig. 29. Sequence C contains the largest 
number of LPVs (cf. bottom panel of Fig. 28), most of them 
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Fig. 29. Log of the second period versus [3.6] diagram for double- 
periodic LPVs. Colors are as in Fig. 22. 
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Fig. 31. Period- amplitude diagram, on a log scale, for mono-periodic 
LPVs according to the sequences of Fig. 25. Colors are as in Fig. 22. 
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Fig. 30. Log of the second period versus [3.6] diagram for double- 
periodic LPVs, according to the sequence in which the first period was 
found. Colors are as in Fig. 22. 



having their first period in sequence D, as noted first by Wood 
et al. (1999). There is also a significant number of stars in se- 
quence C coming from sequences C and E. One can see that 
most LPVs with their first period in sequence C (resp. C) have 
their second period in sequence C (resp. C) (cf. middle panel 
in Fig. 30). This exchange between these two sequences from 
the first to the second period, suggests that they are semi-regular 
variables. Sequence D is one of the most populated, mainly be- 
cause of LPVs with their first period in either sequence A or B. 



There are also stars coming from sequences C and C, in partic- 
ular C-rich and extreme AGB ones. Sequence B, and to a lesser 
extend sequence A, in Fig. 30 are dominated by stars whose first 
period places them in sequence D. There is also a swap between 
sequences A and B from the first to the second period. It is more 
noticeable for semi-regulars with their first period in sequence 
A for which half have their second period in sequence B, than 
for first-period sequence B semi -regulars for which only 18% 
have their second period in sequence A. Sequence E variables 
represent 10% of the double-periodic LPVs. Half of them have 
their first period in sequence D, the other half being distributed 
among the other sequences. 

5.4. Period-amplitude relations 

Figure 31 shows the period-amplitude diagram for mono- 
peiiodic LPVs for each of the sequences defined in Sect. 5.3.1. 
Amplitudes were computed as a peak-to-peak difference in the 
light curves from the R band. From fundamental to high overtone 
mode pulsators, one can see that the mean amplitude decreases 
(cf. Table 6), in accordance with theory, which predicts that fun- 
damental mode pulsators have larger amplitudes than overtones 
ones. Within each sequence, we also observe an increase in am- 
plitude with period depending on the type of LPV as defined 
at the end of Sect. 5.1. The RGB stars have the shortest peri- 
ods and on average the smallest amplitudes, whereas extreme 
AGB stars, when present, have the longest periods and largest 
amplitudes (cf. Table 6). The O-rich and C-rich LPVs have in- 
termediate periods and amplitudes, O-rich LPVs having shorter 
periods than C-rich ones. However, it is impossible to visually 
distinguish these two types in terms of their amplitudes. 

Apart from a few cases, the amplitude for stars from se- 
quences A and B remains below 0.3 mag (log Amp = -0.5). 
As previously mentioned in Sect. 5.3.3, a mix of SRVs and 
OSARGs is possible in these sequences. Sequence C contains 
mainly O-rich and C-rich LPVs with amplitudes between 0.2 
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Table 6. Basic parameters of the sequences of Fig. 23. We assigned one sequence per star according to its first period found. 
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Fig. 32. Period-amplitude diagram, in log scale, for the first period of Fig. 33. Period-amplitude diagram, in log scale, for the second period 



double-periodic LPVs according to the sequences of Fig. 25. 



of double-periodic LPVs according to the sequences of Fig. 25. 



(logAmp - -0.7) and 0.8 mag (logAmp - -0.1), with a few 
extreme AGB stars at amplitudes larger than 1 mag. 

Sequence C contains RGB to extreme AGB stars with most 
of the amplitudes between 0.3 (logAmp = -0.5) and 5 magni- 
tudes (log Am/? - 0.7). 

The RGB stars form the main population of sequence E and 
have a low mean amplitude. No clear linear correlation between 
period and amplitude is visible. However, the lower envelope 
formed by the distribution of RGBs seems to show a decrease 
in amplitude with period up to logf - 2 followed by an in- 
crease in amplitude toward higher periods. At longer periods, 
logf > 2.5, an increase in amplitude is visible for O-rich and 
C-rich LPVs, extreme AGB stars being again among the largest 
amplitude LPVs. 

The shape of sequence D appears very similar to that of se- 
quence E. It could contain two populations: LSPs made of O- 
rich, C-rich, and longer period RGB stars, and stars that indeed 
could be a part of the faint end of sequence E, visible as a group 
of stars with [3.6] > 15 mag and 1.5 < log f < 2 in Fig. 25. 

Figures 32 and 33 show the period-amplitude diagrams for, 
respectively, the first and second periods of double-periodic 
LPVs for each of the sequences defined in Sect. 5.3.1. Overtone 
pulsators from sequences A and B display the same amplitude 
and period properties as for mono-periodic LPVs. 



Most O-rich and C-rich variables from sequence C (exclud- 
ing the extreme AGB stars) have amplitudes smaller than 1 mag- 
nitude, whereas mono-periodic variables extend to larger ampli- 
tudes. In addition they have the same amplitude range as the 
ones in sequence C. This suggests that double-periodic LPVs 
from sequences C and C are semi-regular variables with each 
of their periods in one of these sequences. Extreme AGB stars 
are still among the largest amplitude LPVs and hence are clearly 
distinguished from the other types. Most of them have one of 
their periods on sequence D (cf. Fig. 30), which makes them 
LSP variables. 

Except for RGB stars with the shortest periods, almost all 
double-periodic stars from sequences E and D follow peiiod- 
amplitude relations, the shorter the period the smaller the am- 
plitude. The RGB, O-rich, and C-rich stars have amplitudes 
smaller than one magnitude, and the differences between the 
mean amplitudes of the RGB stars, O-rich, and C-rich AGB stars 
are larger for these two sequences than for the other ones (cf. 
Table 6). Extreme AGB stars are LPVs with the longest periods 
and the largest amplitudes. 



6. Summary 

We have performed a systematic search for LPVs in the variable 
star database from the EROS-2 survey of the LMC. We have 
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shown that applying the criterion for the Abbe value r of r < 0.4 
leads to a selection of smooth light curves over the 6.5 year dura- 
tion of the survey. A cut based on the red giant clump color was 
used to separate our LPV candidates from other types of variable 
stars. 

We have created a catalog of 43 55 1 LPV candidates. Combining 
three period-search methods, we obtain periods for more than 
86% of these candidates. Thanks to data in J, K,, and [3.6] bands 
from 2MASS and Spitzer, we have determined the O-rich or C- 
rich status of our LPV candidates. We retrieved the classical se- 
quences in the period-magnitude diagram studied by Wood et al. 
(1999) and Ita et al. (2004). We have demonstrated that these 
sequences are easier to distinguish using [3.6] data than in the 
near infrared. The main improvement is for the extreme AGB 
stars that are enshrouded into a dust shell whose spectral energy 
distribution peaks around a few microns. 

Our analysis of period-magnitude diagrams has shown that 
mono-periodic LPVs are preferentially distributed on the fun- 
damental mode sequence C and ellipsoidal binaries from se- 
quence E. Double-periodic LPVs have their first period mostly 
distributed in sequence D and their second period in first/second 
overtone sequences C and B. Extreme AGB stars appear as the 
most luminous ones and are mainly found on the fundamental 
mode sequence C of the mono-periodic cases. The C-rich stars 
studied are visible from sequence C up to the highest periods 
in sequence D. Less luminous O-rich AGB stars are distributed 
all along the sequences and are clearly separated from the RGB 
stars by a clear drop in the stellar density at [3.6] = 11. 85 mag. 
The RGB stars are spread over all the sequences and constitute 
the main component of sequence E, which consists of binary sys- 
tems. Our study of period-amplitude diagrams has shown that 
amplitudes tend to increase slowly as we consider longer peri- 
ods from RGB to C-rich AGB stars. Extreme AGB stars have 
the highest amplitudes among LPVs in our sample, well above 
those of the other subtypes. Our analysis of the complete EROS- 
2 variable stars database to extract LPVs remains ongoing work 
that will allow us to perform a homogeneous study of LPVs in 
different metallicity environments. 
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